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Abstract 

We propose, analyze and demonstrate an architecture for scalable cooperative reception. In a cluster 
of TV + 1 receive nodes, one node is designated as the final receiver, and the N other nodes act as 
amplify-and-forward relays which adapt their phases such that the relayed signals add up constructively 
at the designated receiver. This yields received SNR scaling linearly with N, while avoiding the linear 
increase in overhead incurred by a direct approach in which received signals are separately quantized 
and transmitted for centralized processing. By transforming the task of long-distance distributed receive 
beamforming into one of local distributed transmit beamforming, we can leverage a scalable one-bit 
feedback algorithm for phase synchronization. We show that time division between the long-distance 
and local links eliminates the need for explicit frequency synchronization. We provide an analytical 
framework, whose results closely match Monte Carlo simulations, to evaluate the impact of phase noise 
due to relaying delay on the performance of the one-bit feedback algorithm. Experimental results from our 
prototype implementation on software-defined radios demonstrate the expected gains in received signal 
strength despite significant oscillator drift, and are consistent with results from our analytical framework. 
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I. Introduction 

Distributed MIMO (D-MIMO) refers to a broad class of techniques in which a group of cooperating 
nodes emulate a virtual antenna array, in order to obtain performance gains similar to those provided 
by conventional centralized MIMO. While it is difficult to scale centralized arrays to a large number 
of elements due to size and weight considerations (especially at lower carrier frequencies), in principle, 
D-MIMO allows us to synthesize very large apertures using the natural geographic distribution of the 
cooperating nodes, and offer an approach to massive MIMO that sidesteps form factor constraints. The 
opportunistic formation of D-MIMO clusters can also have significant benefits in enhancing range/rate 
tradeoffs, especially in emergency and disaster relief scenarios. However, key difficulties in translating 
D-MIMO from concept to practice are that the cooperating nodes have independent oscillators, each with 
stochastic drift, and that we cannot rely on a regular array geometry in our signal processing algorithms. 
Another important consideration in D-MIMO system design is that we would like our architectures and 
algorithms to scale gracefully as the number of cooperating nodes increases, in order to approach the 
vision of arbitrarily large virtual arrays. In this paper, we address these issues in the context of distributed 
receive (D-Rx) beamforming. 

In a D-Rx beamforming system, a cluster of nodes coherently combine their received signals in order 
to enhance the received signal-to-noise ratio (SNR). In a centralized receive array, depicted in Figure [T(a)1 
this is accomplished by routing signals from different receive antennas along wires, with phase shifts for 
coherent combining applied at RF or IF, or digitally at baseband, after downconversion and analog-to- 
digital conversion. An analogous approach for distributed receive beamforming, shown in Figure |l(b)[ is 
for each node to send its received signal to a centralized processor (typically via a fast local wireless link), 
which then applies the appropriate phase shifts to achieve receive beamforming. With this approach, the 
cooperating nodes do not even have to be synchronized a priori. The centralized processor has access to 
the received signal for each node, and hence can estimate relative frequency and phase offsets and then 
compensate for them when combining the signals. However, this direct approach does not scale to a large 
number of cooperating nodes, since the amount of local communication is proportional to the number of 
nodes. We therefore propose and investigate in this paper an alternative approach which attains scalability 
by using “over the air” coherent combining. 

The proposed architecture (discussed in more detail in the next section) is depicted in Figure |l(c)[ A 
receive cluster of N + 1 nodes wishes to enhance the SNR of a signal arriving from a distant source. One 
of the receive nodes is designated as the final receiver, and the remaining N nodes relay their received 
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(a) Traditional Rx beamforming 



(b) Naive D-Rx beamforming 



(c) Scalable D-Rx beamforming 

Fig. 1. Rx beamforming with traditional MIMO, naive D-MIMO and scalable D-MIMO 


signals, adjusting their phases such that the relayed signals combine coherently at the designated receiver. 
This converts the task of D-Rx beamforming on the “long link” from source to receive cluster into one 
of distributed transmit (D-Tx) beamforming on the “short link” between the relays and the receiver. 
This allows us to leverage prior work on scalable D-Tx beamforming, while adapting to features and 
impairments peculiar to our relay-based architecture, in order to attain a scalable D-Rx system. 

The key features and contributions of this paper are summarized as follows: 

1) Scalable architecture: By using amplify-forward relays in the receive cluster, we ensure that local 
communication overhead does not blow up with the number of cooperating nodes. We use a provably 
convergent one-bit feedback algorithm for distributed transmit beamforming in order to ensure that 
the relayed signals add up coherently at the receiver. 

2) Implicit frequency synchronization: For time division between the long link and the short link as 
considered here, the frequency offsets of the relays “cancel out” on the long and short links, hence 
there is no need to synchronize the relays in frequency. 

3) Analytical performance characterization: While frequency synchronization is not required, the 
delay between message reception on the long link and message relay on the short link leads to 
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phase errors accumulating because of frequency and phase drift. We characterize the statistics of 
such phase errors as a function of local oscillator (LO) parameters, and then provide an analytical 
framework for determining their effect on the one-bit feedback algorithm. Our analysis matches 
closely with Monte Carlo simulations, and shows that when the phase error gets large, only a 
fraction of the expected beamforming gain is achieved. 

4) Proof of concept: The proposed architecture is implemented on a software-defined radio testbed, 
showing that the expected gains can be achieved with up to four relay nodes. By relaxing the system 
parameters, we are able to observe the performance of the one-bit feedback algorithm under larger 
phase errors experimentally, thereby verifying the insights from our analytical framework. 

Related work: Many information-theoretic analyses, ranging from three decades back [1, 2] to the 
present [3-5], rely on the concept of cooperative wireless communication, without explicitly addressing 
the fundamental bottleneck of synchronization between cooperating nodes. However, there has been 
significant recent progress on the problem of distributed synchronization, most of it in the context of 
distributed transmit beamforming. Closed-loop synchronization techniques include explicit feedback 


one-bit aggregate feedback [8^100, irn 


21-< 

me- 

12, 

13] 


1311 or two-way synchronization [ 1 


imp 

li 


icit feedback using reciprocity 111 lh . round-trip synchronization 


15fl. Each of these synchronization techniques offers different 


trade-offs between complexity, coordination overhead and scalability to larger networks (see the review 
paper in |16] and the discussion in [17]). Distributed MIMO techniques have also been investigated in 
the context of “coordinated multipoint” (CoMP) capabilities for 4G-LTE systems, where multiple base 
stations act as a distributed antenna array 118], Recent work has shown that the overhead for D-Rx 
beamforming using the architecture in Figure [T(b)1 can be reduced by heavily quantizing the information 
exchanged. However, the local communication overhead still scales up with the number of cooperating 
nodes, unlike our proposed approach. 

There is also by now a significant body of research in analysis and simulation of amplify-and- 


forward relaying: in [2.QD, the receiver broadcasts a single bit of information to each relay indicating 


whether it should participate in the communication, thereby selecting a set of relays who happen to be 


combining quasi-coherently; [21] considers network beamforming where each node has perfect channel 


state information; [22] considers network cooperation where each node computes beamforming weights 


based on local channel information only; and [23j, 


24] 


scheme based on partial channel state information. In [190, nodes forward their message over a local area 


proposes a robust collaborative beamforming 


network, and it is shown that by selecting a subset of the receiving nodes, the performances are close to 
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the ones obtained with optimal receive beamforming. 

While beamforming is implicit in the concept of ideal amplify-forward relay, other than our own 
prior conference paper 12511 . to the best of our knowledge, the present paper is the first to identify the 
critical importance of amplify-forward relay in providing a scalable architecture for distributed receive 
beamforming, to model in detail the synchronization issues in implementing this approach, and to provide 


a proof of concept via our testbed. A significant contribution of this paper beyond [25] is that we provide 


a detailed analytical framework for the impact of phase errors on the one-bit feedback algorithm. We 


also extend the prototype in 12511 to a larger receive cluster, and provide a more detailed characterization 


of the effects of the system parameters on experimental performance. 

Distributed MIMO techniques have also been demonstrated with a wide variety of experimental 


prototypes. In 


I sl Q 


2711 . distributed transmit beamforming prototypes using wired feedback channels 


were presented for both RF and millimeter wave frequencies. A first fully wireless setup was presented in 


[28|], but this setup still used analog signaling for the feedback channel. A D-Tx beamforming prototype 
using full digital signaling and an extended Kalman filter for frequency synchronization was presented 


in [29H32]. While the preceding papers focus on distributed transmit beamforming, multiuser distributed 


MIMO has been demonstrated in [33]. However, the latter uses dedicated wired backhaul links to distribute 
information and feedback throughout the network. Similarly, CoMP experiments for 4G cellular systems 
have to date relied on dedicated backhaul links with low latency and high throughput, as well as assuming 
uninterrupted GPS connections to synchronize the various base stations. To the best of our knowledge, 
the present paper is the first work to analyze and prototype all-wireless distributed receive beamforming. 

Outline of the paper: Section [TT] discusses our system model, the challenges of achieving synchroniza¬ 
tion at multiple levels in such a system, and points out the implicit frequency synchronization achieved 
by the design considered here. In Section [nil we characterize the phase error accumulating due to 
relaying delay. Section [[V] presents the phase synchronization algorithm and studies the performance of 
the this algorithm under the presence of phase error. Finally, Section [V] presents the implementation of 
our architecture on an experimental testbed, and shows some of the results obtained with our prototype. 
The prototype is run with different parameters to confirm the theoretical insights of the previous sections. 
Section [VI] contains our conclusions. 


II. System Model 

While there are many possible design choices for the scalable D-Rx architecture depicted in Figure [T(c)[ 
the specific choices in our modeling, analysis and prototyping are as follows. We focus on narrowband 
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signaling, with channels modeled as complex gains. Each relay node receives the signal from the distant 
transmitter over the long link, applies a phase shift to the received signal, and forwards it to a central 
receiver over a short link. All relays forward the message to the receiver over the same frequency band, 
with over-the-air combining at the receiver. We assume that the same frequency band is used, via time 
division, for the long link and the short link, which enables implicit frequency synchronization. The 
relays adapt their phases using the one-bit feedback algorithm 11011 . which also formed the basis for our 
prior prototyping of distributed transmit beamforming in 13ill . Note that the relays can be implemented 
in RF or in baseband, assuming sufficient ADC resolution. 

The key challenge in coherent combining at the receiver is that signals emitted from relays with 
independent clocks and oscillators must line up. Three levels of synchronization are required: frequency, 
phase and timing. In the following, we describe our approach to each of these problems. 


A. Implicit frequency synchronization 

In D-MIMO systems, each terminal derives its RF signal from its own local oscillator (FO) which 
carries a small but non-zero frequency offset with respect to those of the other nodes. This results in 
FO frequency offsets that can range from a few Hz to tens of kHz for poor quality oscillators. If the 
signals forwarded by the relay nodes to the central receiver have significant frequency offsets, the total 
signal at the receiver will exhibit constructive and destructive interference patterns. In order to avoid 
such behavior, it is important that the signals arriving from the different relay nodes have no frequency 
offset with respect to one another. Fortunately, we get this for free in our architecture, as discussed in 
the following. 

While we have transformed our D-Rx problem into one of distributed transmit beamforming along 
the short link, a key difference from pure transmit beamforming context as in [31] is that we obtain 
implicit frequency synchronization by virtue of time division between the long and short links. To see 
this, suppose that the FO frequency offset between the transmitter and final receiver is denoted by fo, and 
that the offset between the transmitter and the i- th relay is given by fn. The FO frequency offset between 
the i-th relay and the final receiver is then given by fa - fo - fa. Because each relay node forwards the 
signal from the receiver without decoding, the forwarded signal still contains the FO frequency offset 
of the first link. It then follows that all messages arriving at the receiver will carry the same total FO 
frequency offset of fn + fa = fo, which is independent of the relay node i. Of course, in practice, FO 
phase and frequency drift causes an accumulation of phase error over the time between reception on the 
long link and relay on the short link. For the low-quality oscillators in our software-defined radio testbed, 
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such effects are significant. In the next two sections, therefore, we model such phase errors and analyze 
their impact on performance. 


[10], 


B. Phase synchronization with the one-bit feedback algorithm 

Phase synchronization is achieved by using the well-known one-bit feedback algorithm llOtl . which is 
an iterative, scalable technique that has been shown to converge to an optimum both theoretically and 
experimentally. The algorithm is a stochastic hill climbing procedure, with one iteration per cycle. Each 
relay adds a random phase perturbation to its current phase, and forwards its received message. The 
receiver monitors the received signal strength (RSS) of the received, over-the-air combined message, and 
broadcasts back a single bit to the relays indicating whether or not the RSS has increased (compared to 
the earlier maximum RSS). If the RSS has increased, each relay keeps the last random phase perturbation 
and the new maximum RSS becomes the current RSS; if not, each relay discards its previous random 
phase perturbation and the maximum RSS remains unchanged. In the next cycle, the whole process is 


repeated. It has been proven theoretically and with numerous prototypes 
algorithm, the RSS quickly converges to its maximum value. 


fl 


28 


3111 that with this 


This ideal version of the one-bit feedback algorithm does however suffer from LO phase drift, which, 
unfortunately, occurs in any real world implementation. Consider the following scenario: at some point 
during the one-bit feedback algorithm, a combination of phases at the relay nodes is achieved that 
maximizes the RSS. Subsequent iterations of the algorithm will therefore not be able to improve on the 
previously obtained maximum. However, phase noise will cause the relays’ phase offsets to drift away 
from their ideal values. Over time, the phase offsets will accumulate more random error, the RSS will vary 
randomly at lower values, and the one-bit feedback algorithm will fail entirely. To avoid this undesirable 
behavior, the following change is applied to the one-bit feedback algorithm. Instead of comparing the 
RSS in a given cycle to the maximum RSS achieved in a all previous cycles, the current RSS is compared 
to the maximum RSS in the K previous cycles (where K is a finite number). The algorithm will then 
be able to recover from phase drifts, as the RSS at each cycle is only compared to the K previous RSS 
values. If at some point a combination of phases is obtained with an RSS that cannot be outperformed, 
K cycles later this maximum RSS will be removed from the past RSS memory, making the algorithm 
robust against phase drifts. 
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C. Implicit timing synchronization 

The third type of synchronization is temporal: the packets forwarded by the relays must arrive simul¬ 
taneously at the receiver to avoid the effects of inter-symbol interference (ISI). In this work, we assume 
a narrowband channel, so that the differences in propagation delays are significantly smaller than the 
inverse bandwidth and can be neglected. In that case, determining packet boundaries on the long link 
accurately and delaying the forwarded packet for a fixed amount of time is sufficient to obtain precise 
message alignment at the receiver. Thus, we use the timings of the messages received by the relays 
on the long link to provide implicit timing for the messages relayed on the short link. For wideband 
dispersive channels, more sophisticated strategies are required to handle ISI (e.g., OFDM with frame 
synchronization across relays), but this is left as a topic for future work. 


III. Modeling phase errors 

We now model the phase error accumulated due to relaying delay. Suppose the transmitter sends a 
message at time instant t\. We know that the LO frequency offset satisfies the following condition: 

fn(h) + fi2(h) = fji(h) + f j 2 (ti) = /o(fi) Vz, j (1) 

Flowever, in reality there is some time delay T d before the relays forward the message to the final receiver. 
The LO frequency offsets of the different relay nodes will drift independently, such that £Q| is no longer 
satisfied, which can be expressed by 

fn(h) + /t 2 (f 1 + T d ) * + fj^ti + T d ) \/i,j 


As a result, the signals arriving at the receiver from different relays now have different LO frequency 
offsets, leading to fading in the total message. 

We now model the effect of LO drift in more detail. We show that it causes a zero mean Gaussian 
phase error, and compute the variance of the phase error as a function of system parameters. 

LO model: We begin by reviewing the following state space LO model 04 35]. Suppose (P(t) and 
io(t) are the LO phase and (angular) frequency offset at time t. The LO phase and frequency offset at 
time t + T are given by: 


Ht + T) 

oj(t + T') 


1 T 
0 1 


4Kt ) 

u(t) 


+ n(T) 


where n(T) is the process noise vector, distributed as n ~ jV(0, Q(T)), which causes the LO phase 
and frequency offset to drift from their nominal values. The process noise covariance matrix Q(T) is 
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modeled as [34 


35|] 


Q (T) = u 2 A 


T 

0 

, 2 2 

rj- 13 rp 2 

T T 

0 

0 

+ uj cQ2 

T 2 rp 

.2 1 


where uj c is the carrier frequency and the parameters q\ and g 2 are the process noise parameters that 


coiTespond to white frequency noise and random walk frequency noise, respectively. 


Intra-cycle drift: Consider the signal at the final receiver shown in Figure [2] At time t\, the distant 
transmitter sends its message to the relays (received with low amplitude), and at time f 2 , the relays 
amplify and forward the message to the receiver (received with higher amplitude). The phase of the 
relayed signal at the receiver is then given by 


4>i(t 2 ) - + ^* 2 (^ 2 ) 


^2 

= 4>ii(ti) + <t>i2(ti) + T d LUi 2 (ti) + n (j)i2 (T d ) + J n Ui2 (T)dr 

tl 


where <j>n and 0 l2 denote the phase offset of the long and the short link, 0i denotes the phase offset of 
the global link (from transmitter to receiver) through relay i, is the LO phase noise over the short 
link and n Wi2 is the LO frequency noise over the short link. The variance due to phase noise is then given 

by ujcqfTd, whereas the variance due to the LO frequency noise is given by cu c q 2 -y. The latter follows 

b 

from the well-known result that for a standard Wiener process W, we have that J Wtdt ~ jV(0, 

a 


The total phase error variance due to intra-cycle drift is thus modeled as a zero mean Gaussian variable 
with variance 


rp3 

2 2 2rp , 2 2 ± d 

a <P - ^c/h Td + w c q 2 — 


( 2 ) 


The increase in variance with relaying delay Td is intuitively reasonable: the longer the relays wait before 
forwarding the message, the larger the phase errors accumulated due to the independent LO drifts at the 
different relays. 
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Fig. 2. Received signal over two cycles. 
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Inter-cycle drift: Figure [2] shows several relay cycles. At time t\, the transmitter sends a first message, 
which is forwarded by the relays at time £ 2 - The receiver then computes its feedback message that is 
returned to the relay nodes; essentially, this message controls the phase 7 ,; that is applied by each relay, 
such that the total phase </>ii(fi) + 0 , 2 ^ 2 ) + 7 i is equal for all links. Now consider a second cycle of the 
system: at time 1 3 , the transmitter sends another message, which is forwarded by the relays at time £ 4 . 
The total phase of the second message for link i is given by 

4>i(U) = 4>n(t 3 ) + 4>i2(t4,) +7i 

u 

= <^1(^3) + <t>ii(h) + T d u i 2 {t 3 ) +n (j> . 2 (T d ) + f n Ui 2 (T)dr +7; 

£3 

u 

= +(j)i2(t 3 ) +T d {u} i 2 {t 1 ) +n UJi 2 (T c )) + n < f >a (T d ) + J n Ui 2 (T)dT +7; ( 3 ) 

^3 

where T c is the cycle time. It can be seen that between t 1 and t 3 the LO frequency offset 072 has drifted, 
potentially causing the relay phase shift 7 * to become outdated. Obviously, if the cycle period gets longer, 
this drift will be more severe, resulting in a larger phase error. From ©, the total phase error due to 
both intra- and inter-cycle drift can be characterized by a zero-mean Gaussian process with variance 

rnS 

2 2 2/t - ! 2 2 d , 2 2r7 _ i2r7~i ^a\ 

- u c qiTd + w c q 2 — + (4) 

The first two terms of the equation are the same as for © and represent the intra-cycle drift, whereas 
the last term represent the inter-cycle drift. It can be concluded from © that, for fixed LO parameters 
q\ and g|, both the relay delay time T ( j and the cycle time T c should be kept low in order to maintain 
small phase errors. 

In the next section, we analyze the effect of phase errors on the performance of the one-bit feedback 
algorithm. 


IV. One-bit feedback with phase errors 

In Section [nil we showed that implicit frequency synchronization results in a zero-mean Gaussian 
phase error for each relay. Even for poor quality LOs, these phase errors are well below 27 t, and do not 
affect frequency synchronization. However, they do affect phase synchronization, and in this section, we 
develop an analytical framework to characterize its impact on the one-bit feedback algorithm (with RSS 
memory of length K\ see description in Section lll-BI ). In Section |V| the effects of increased phase errors 
is investigated experimentally. 
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The key steps in our derivation are the following: 

1) We determine the probability that the noisy phases cause an RSS increase, conditioned on the fact 
that the noiseless phases cause an RSS increase, over a single cycle of the algorithm (Section [IV-AI ). 

2) In the case of finite memory, we define K states Sk as being the situation where the maximum RSS 
was achieved k cycles ago, and model the transitions between states as a Markov chain. Following 
arguments similar to the ones in Section IIV-A1 we compute the phase transition probabilities and 
derive the probability of being in each state Sk (Section IlY-BI ). 

3) We compute the RSS drift - the expected deviation in RSS at each iteration in the algorithm - when 
operating with noisy phases, both for positive and negative feedback, conditioned on the current 
state (Section IlY-BI ). 

4) Using the probabilities of being in each state Sk and the RSS drift conditioned on state Sk, we 
determine the total RSS drift when running the one-bit feedback algorithm with noisy phases 
(Section IIV-B b . 


A. One-bit feedback algorithm with Gaussian phase noise over a single cycle 

We start by investigating the effect of Gaussian phase noise over a single cycle. For N nodes, we 
define the following normalized RSS values 


V = 


ys = 


N 

El 

2=1 




1 N 

_ y | ^ (<£<+*<) | 
1 


Vn = 


2=1 

N 

El 


VSn 


1 ' i =1 

1 N 

-Y 

Nti 




Jifa+Si+n'A 


that represent the noiseless RSS before random phase perturbation, the noiseless RSS after random phase 
perturbation, the noisy RSS before random phase perturbation and the noisy RSS after random phase 
perturbation, respectively. In these equations, <5* is the random phase perturbation applied at each node, 
and m, n\ are the phase noises (before and after random phase perturbation) at node i. 

The problem at hand is then to determine when the one-bit feedback performs as well in the noisy 
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case as in the noiseless case. To begin, we define 

U ~ V5n — Vn 

v = ys - y 


which denote the RSS increment after a phase perturbation in the noisy and noiseless cases, respectively. 
In addition, let us define the state of the one bit feedback process as the value y of the RSS in the 
noiseless case. Conditioned on the feedback process state, the probability of successful detection of a 
phase improvement/deterioration is then 


P[f7 > 0|V > 0,2/] 
F[U < 0|V < 0 ,y] 


Clearly, the random variables U and V are not independent, so to determine IP [U > ()\V > 0, y] and 
P [U < 0\V < 0 . y], one must first determine the joint probability distribution of (U,V). To proceed, we 
leverage the results derived in Conjecture 1 in 11 1 Oi l, where statistical mechanics arguments were used to 
determine that that the distribution of the phases around their mean follow an Exp-Cosine distribution. 
Given this framework, we show in Appendix lAl that the joint probability distribution of (U. V) can be 
approximated by the following bivariate Gaussian distribution 


u 


/ 

(XSn - Xn)y 


t-xL -PsnK(y)+l-xl-PnK(y) 

Xn-XSrvXS~p2SnK'(y) 





2N 

2N 


V 


1 

5 * 

T—1 

1 

'O 

>< 

_1 

5 

Xn-X6nXS~P2S n K(y) 

i-xs-pf> K (y) 



\ 


L 2 N 

2 N j 



where X5 = E[cos( 6 j)], Xn = E[cos(rc,)] and xSn = E[cos(5; + n')], where p s = xj ~ E[cos(25;)], p n = 

X 2 -E[cos(2n;)] and p Sn = x 2 n -E[cos(2(5j+ nj))], where p 2 s n = XSnXS -E[cos(2^ + «*)], and where 
N 

the term n(y) - T, E[e 2j ®' ] depends on y only, and can be approximated by n(y) = eT 4i 1 ^ y ‘ > for large 
i=1 

y. From ©, one can easily compute the probability that the one-bit feedback algorithm is successfully 
able to detect a phase improvement P [ys n > Vnlus > !)■, y] or a phase deterioration P [yg n < y n \ys < y, y] 
in the noisy case. 


The joint distribution © has been compared with Monte-Carlo simulations and is shown to match the 
simulations very well for as few as 10 nodes. In these simulations, the random phase perturbation was 
chosen uniformly from the discrete set {-10°,+10°}, as to match our experiments in Section [V] Note 
that other distributions for the random phase perturbations yield similar conclusions. The phase noise is 
drawn from a zero-mean Gaussian distribution with variance a\. It can be seen in Figure [3] that, for low 
phase noise, U and V are highly correlated: a successful decision in the noiseless case gives way to an 
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identical decision in the noisy case. As the phase noise increases, U and V become less correlated, and 
the probability that the one-bit feedback algorithm makes a correct decision (with respect to the noiseless 
case) becomes lower. 


N = 10 nodes — y = 0.8 — 5. ~ [-10°, 101 



0.2 


( c > a n = 10° (d) cr n = 15° 



- 0.2 - 0.1 0 0.1 0.2 - 0.2 - 0.1 0 0.1 0.2 
Vsn-Vn V^n 


Fig. 3. Comparison between Monte-Carlo simulations and theoretical model for various phase noise variances and N = 10 
nodes. The random phase perturbations are chosen uniformly from the discrete set {-10°,+10°}. 


B. One-bit feedback algorithm with Gaussian phase noise and an RSS memory of length K 

As mentioned, a practical implementation of the one-bit feedback algorithm requires the use an RSS 
memory of finite size. We therefore consider the case where the one-bit feedback algorithm compares 
the current RSS at time t with the RSS of the past K cycles. In the following framework, we define the 
state Sk as the state where the maximum RSS during the I\ previous cycles was obtained during cycle 
t - k, as shown in Figure [4j At each time instant, we are in one of the K possible states S^. 

We are interested in computing the total RSS drift, that is the average RSS increment at time t 
conditioned on the current RSS level y. We make the following simplifying assumption: if we are in 
state Sk at time t, the RSS drift is statistically independent of the feedback before time t - k. In other 
words, the RSS drift is only dependent on the feedback obtained between time t - k and time t. This 
can be verified as follows: imagine the following time instants tq <T\ <t-k < T2 <t, where tq and t-k 
correspond to two time instants when there was a phase update. Instant tq is the reference point at time 
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Fig. 4. Different states of the system for a memory of size K = 4. 


Ti, and t - k is the reference point at time 12 and time t. Since the random phase perturbations and phase 
noises are independent across iterations, we can show that the covariance between U Tl and U T2 and the 
covariance between Vt and U Tl is equal to (using the notations defined in Appendix lAl ) 


Cov [U Tl , Ur 2 ] — IE Tl %9\,n,TQ S ) Sn,r 2 j 0 

CoV [Vj, C^r 2 ] — IE ~ X9\,n,Ta) j — 0 


which shows that the RSS drift at time t is independent of the feedback before t - k. 

Given the previous assumption, the (noiseless) RSS drift conditioned on the current state ,S'/. is given 
by 

Drift(RSSt\Sk,y) = IE [V(|feedback since t - k, y] 

The total RSS drift with a memory of length K is then given by: 

I< 

Vnft(RSS t \K, y) = £ P (S k \y) ■ Dnft(RSS t \S k ,y) 

k =1 
K 

= Y,P(S k \y) ■ [P (U t > 0\S k )Dnft(RSS t \S k ,U t > 0,y) + F(U t < 0\S k )Dnft(RSS t \S k ,U t < 0,y)] 

k =1 

( 6 ) 

We begin by determining the probabilities P(S’/ c |y) of being in a given state k. We model the state 
transitions with the Markov chain shown in Figure [5J In each state, a positive feedback will bring the 
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Fig. 5. Markov chain modeling the state transitions for a memory of size K = 4. 


algorithm back in state S\, whereas a negative feedback will cause the algorithm to transition from state 
S/. to state Sfc+i- For the final state Sk, a negative feedback will cause to algorithm to transition to any 
state with uniform probability, as the noisy RSS between sates t, t - 1 ,t - K + 1 are i.i.d. The Markov 
chain state transition matrix is given by (for a memory of length 4) 


Pl 

1 ~P\ 

0 

0 

P2 

0 

1 ~P2 

0 

P3 

0 

0 

1 - P'i 

P4 + -r 1 

l-Pi 

4 

1 Pi 

4 

1 — jP4 

4 


By looking at the stationary distribution n of the Markov chain, we can determine the probability of being 
in each state. The stationary distribution ir is the left eigenvector of P that correspond to the eigenvalue 


A = 1: 


7T T P = 7 T T 


( 8 ) 


which will depend on the state transition probabilities pf,. and the memory length I\. 

We now determine the state transition probabilities pk, which define the probability of a positive 
feedback at a given state Sk- The probability /y. is defined as 

p k =P [U t > 0|S fc ;y] 

= p [Ut> 0\U t -i < 0,..., Ut-k+i < 0; y] (9) 

Using similar arguments as in Section HV-Al we model the random variable \Ut, Ut-i, ■■■, Ut-k+i] as a 
multivariate Gaussian random variable. The mean and variances of each term have already been computed 
in Section ITV-Al and the covariance puu between two variables Up and Up is given by (using arguments 
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similar to the one in Appendix |0 

PUU — ^ ~ l ) (®9t ,Sn,t 2 ~~ -P9i,<5n,t-z)] 

= ^ \ X 9K,5n,t-l\ 

= 1 ~ XSn ~ PSnKjy) 

2N 

Note that this term is identical for all possible values of t\ and t- 2 - Once we know the mean and 
full covariance matrix of the multivariate Gaussian random variable [Ut, Ut-i,U t _k + i] , we can use 
determine the probability p given in ([h]) by computing the multivariate cumulative density function. By 
tilling the values of in (0 and one can compute the probability P(S'/-) of being in state k. Figure [6] 
shows the probability of being in each state for N = 10 nodes and a memory of size K = 4, and compares 
our theoretical model with Monte-Carlo simulations. It can be seen that there is a good match between 
theory and simulations. 


N = 10, K=4, o n = 1” 



Fig. 6. Probability of being in each state, for N = 10 nodes, a memory of size K = 4, a random phase perturbation of 
(5 ~ W[-10°, 10°], and a phase noise of <j n = 1°. 


We now compute the RSS drift terms of equation ©. The drift term for positive feedback is defined 
as 

Dnft(RSS t \S k , U t >0;y)=E [V t \U t > 0, t/ t _i < 0,..., U M < 0; y] (10) 

We again consider the Gaussian multivariate distribution \Vt, Ut , Ut- i, •••, Ut-k+ 1 ]■ The mean and variances 
of all the terms have been computed previously, as have the covariance between and Ut 2 and the 
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covariance pjjy between Ut and V*. The covariance pu T v t between Vt and U T is given by (for t-l < r < t) 

PU T ,V t ~ ® — 3'9t i 6n 1 t—l) 

= IE — •X'9\,5n,t-l )]E[x SHj< 5 jt ] 

= 0 

The covariance matrix of the multivariate Gaussian distribution is then given by 



a v 

Puv 

0 

... 0 


Puv 

a u 

Puu 

••• PUU 

s = 

0 

Puu 

a u 

••• PUU 


0 

Puu 

Puu 

- a u 


Computing the RSS drift term in (flOl) then boils down to determining the mean of a truncated multivariate 
Gaussian distribution, which can be done efficiently with Monte-Carlo integral computation. 

For the case of negative feedback, and for all states k = 1 ,K - 1, the RSS drift is equal to zero, as 
the one-bit feedback algorithm does not change the combination of phases: 

Drift (RSS t \S k , U t < 0; y) = E [V t \U t < 0, U t . i < 0, U t . k+1 < 0; y] = 0 (11) 

Note that this equality does not hold for k = K. If Ut < 0, the previous maximum RSS then “slips out” of 
the past RSS memory, and in the next iteration a new maximum RSS will be considered. By definition, 
this will cause a change in true RSS drift, which can be expressed as 

Dnft(RSS t \S K , U t <0;y) = ¥ [Max. noisy RSS at t - r] E [V t - T \U t - T < 0 , {U t -i < U t - T } K -i<T<t,l*^ v] 

= E [V t \U t < 0, U t - 1 < U t ,..., U t - K+ 1 < U t -y] (12) 

where this last equation is obtained by symmetry arguments. Equation (fl2l) can be determined by yet 
another multivariate Gaussian distribution with variates {Vt, Ut, Ut-i,Ut-K+i] 1 ■ The elements of the 
mean and covariance matrix of this multivariate Gaussian distributions have all been computed previously, 
and the solution of (IT2T) can be obtained through Monte-Carlo integral computation. 

By combining the different RSS drift terms (flOl) . (fTTb and (fl2l) into equation ©, the total RSS drift 
can be computed. Figure [7] shows the RSS drift for various phase noise values, both with Monte-Carlo 
simulations and with our theoretical model. It can be seen that there is a good correspondence between the 
simulated and the theoretical curves. It should be noted here that to obtain stable Monte-Carlo curves, a 
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large number of realizations must be generated, making the computational requirements quite expensive. 
In contrast, our theoretical analysis is fairly efficient and does not require significant computation power. 
From Figure [7] it can be seen that when the phase noise becomes larger than the random phase perturbation, 





Fig. 7. Theoretical and simulated RSS drift for N = 10 nodes, a memory of size K = 4 and a random phase perturbation of 
5 ~ U\- 10°, 10°], for various phase noise values. 


the RSS drift eventually becomes negative. This means that in steady-state, the RSS will converge to a 
value of y smaller than 1, and not achieve the maximum possible RSS. For a phase noise of a n = 15°, 
the RSS will only reach 80% of the maximum achievable RSS. 

Our analysis is confirmed by the simulation results in Figure [8] Here, the normalized RSS is plotted 
versus time when running the one-bit feedback algorithm. The normalized RSS has been averaged over 
100 simulation runs, and the first 1000 iterations are not plotted in Figure [8] (in order to focus only on 
the steady-state convergence values). It can be seen that the normalized RSS does converge at a value 
that is predicted by the zero-crossing of the RSS drift in Figure [7] In Section [V] we will show that our 
experimental testbed starts failing when the phase noise gets larger than the random phase perturbations. 


V. Experimental demonstration 
A. Software-defined radio testbed 

The proposed architecture was implemented on a software-defined radio testbed using six USRP RF 


and baseband boards [36]. We use a mix of USRP-2 and USRP-N200 baseband boards, and WBX 50- 
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Fig. 8. Normalized RSS when running the one-bit feedback algorithm, averaged over 100 runs. It can be seen that as the phase 
noise increases, the steady-state average RSS decreases. The average normalized RSS can be predicted by the zero-crossing of 
the RSS drift in Figure [7] 


2200 MHz RF daughterboards. Each USRP was connected to a host laptop that performed the computation 
using GNU Radio software I 37 I 1 . Our software is available for download online B. 

One USRP was used as a transmit node (sending packets that contain only a pilot tone), one USRP was 
used as the final receiver, and up to four USRPs were used as relay nodes. A block diagram of the relay 
nodes is shown in Figure |9(a)| The relay nodes receive the message from the transmitter, add a phase 
shift to the received message, and wait for a fixed amount of time T,j before amplifying and forwarding 
the message to the final receiver (this delay 7)/ needs to be identical for all relays for the messages to 
add up without ISI). The message is forwarded over the same frequency band as it is the one used by 
the transmitter, which was set to 908 MHz. Additionally, the relay nodes also listen to the feedback 
message over the feedback channel (at 928 MHz) to determine the phase shift to be applied using the 
one-bit feedback algorithm. The sample rate of all the nodes was set to 200 kHz. The receiver, shown 
in Figure |9(b)} receives the combined messages of all relay nodes, computes its single bit of feedback 
and broadcasts this single bit back to all relays over the feedback channel. The single bit is embedded 
in a GMSK-modulated packet. 

B. Experimental results 

In this section we present results obtained with our experimental prototype. The prototype was run 
in an indoor environment, with a distance between transmitter and relay/receiver node of approximately 
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(a) Relay node 


Rx Tx 

antenna antenna 



(b) Receive node 

Fig. 9. Block-diagram for relay and receiver node. 

5 m. The nodes were kept static during the experiments, and there was little movement around the testbed 
to limit the effects of dynamic fading. The random phase perturbation are chosen randomly from the 
discrete set {-10°, 10°}. Figure ITOl shows the received signal during a single cycle for two relay nodes. 
The relay delay time was set at 7j/ = 10 ms, and the cycle period at T c = 50 ms. In subligure (a), no 
relay is activated: only the message from the (distant) transmitter is observed, with low amplitude. In 
subfigures (b) and (c), only relay 1 or relay 2 is activated. After the original message from the transmitter, 
the (stronger) message from the relay can be observed. Finally, when both relays have been activated 
and convergence of the one-bit feedback algorithm has been achieved, it can be seen in subfigure (d) 
that the relayed packets from both relays add up coherently. The amplitude of the relayed packets is then 
equal to the sum of the amplitudes of the individual relayed packets. Note that in all figures there is a 
noisy signal after the relayed packets. This corresponds to the self-interference created by the receiver’s 
feedback message to the relays, in an adjacent frequency band, and can be ignored. 

Figures [TT] and [12] show the mean amplitude of the relayed packets only, over longer amounts of time, 
using 3 relays and 4 relays, respectively. It can be seen that the amplitude of the combined relayed 
messages correspond to the sum of the amplitudes of the individual relayed messages. Also, it can be 
observed that, once the one-bit feedback achieves convergence, the amplitude of the relayed messages is 
stable at its maximum value. Thus, the phase errors due to LO drift are being successfully handled by 
the one-bit feedback algorithm. In Figure [121 the steady increase in RSS can be observed when the 4th 
relay is turned on. A few iterations were necessary for the RSS to converge to its maximum value. It 
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Fig. 10. Received signal during one cycle of the setup with (a) no relays, (b) relay 1, (c) relay 2 and (d) relays 1 and 2 
activated. 



Fig. 11. Mean amplitude of the relayed packets with 3 relay nodes. 


can also be seen in Figures ITTHT21 that there are slight dips once the RSS has converged to its maximum 
value. This is because the one-bit feedback algorithm continues running even after the RSS has achieved 
its maximum value, causing the phases to misalign and realign over time. An easy improvement would 
be to reduce the size of the random phase perturbation applied at the relays once the RSS converges to 
its maximum value. 
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Fig. 12. Mean amplitude of the relayed packets with 4 relay nodes. 


In Section [Hi] it was determined that increasing the relay delay time Tj and the cycle period T c would 
result in an increasing phase error. In addition, it was argued in Section [IV] that if the phase error becomes 
large (with respect to the size of the random phase perturbation), the RSS drift will decrease and the 
one-bit feedback algorithm will not be able to maintain the amplitude at its theoretical maximum. To 
verify these predictions, our experimental testbed was run with different values of T c and T ( j, as shown 
in Figure [13] The setup was run with two relay nodes, and the random phase perturbation on both relay 


nodes was 10°. The LO parameters of our testbed were estimated previously [3111 as q{ = 8.47 x 1CT 22 
and (J 2 = 5.51 x 10 -18 . The corresponding phase error standard deviation, computed using ([4]). is given in 
the title of the subfigures. For each test, we first waited for a period of time long enough that the one-bit 
feedback algorithm could be expected to converge. The red line represents the (normalized) maximum 
possible RSS (based on the measured amplitudes of the relayed packets when the relays are turned on 
individually), and the blue line corresponds to the (normalized) measured RSS of the relayed packets 
when both relays are turned on, after convergence of the one-bit feedback algorithm. It can be seen that 
once the phase error standard deviation becomes significant with respect to that of the random phase 
perturbation, the one-bit feedback algorithm has trouble converging, and the RSS has trouble maintaining 
its maximum value. 
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Tc = 50 ms — Td = 20 ms — a n = 3.9° Tc = 100 ms — Td = 40 ms — a n = 10.5° 

1 I -.. - - 1 I— < V » 

'0.5 ^ 0.5 

0 

0 5 10 0 5 10 

Time (s) Time (s) 

Tc = 50 ms — Td = 25 ms — a = 4.9° Tc = 100 ms — Td = 50 ms — a = 13.3° 




1 

' 0.5 
0 

0 5 10 0 5 10 

Time (s) Time (s) 

Tc = 50 ms — Td = 30 ms — a = 5.9° Tc = 100 ms — Td = 60 ms —a = 16.1 ° 


1 Ty*v^» 

'0.5 

0 - 


0 5 10 

Time (s) 


; . f 


0 5 10 

Time (s) 


Fig. 13. Mean amplitude of the relayed packets when varying Td and T c . 


VI. Conclusions 

Starting from the observation that over-the-air combining using amplify-forward relaying provides a 
scalable approach to distributed receive beamforming, we have proposed an architecture for achieving 
the synchronization required for the relayed signals to cohere at the receiver. An attractive feature of the 
time division (between long and short links) approach considered here is that frequency synchronization 
comes for free. We have demonstrated this architecture using a software-defined radio testbed, and report 
experimental results achieving the receive beamforming gains predicted by theory. We also model and 
analyze the potential performance degradation due to phase errors accumulating due to LO drift. We 
provide an analytical framework, verified via Monte Carlo simulations, which estimates the degradation of 
the RSS attained by the one-bit feedback algorithm with finite memory in the presence of phase errors. A 
key insight, also verified experimentally, is that significant performance degradation occurs if the variance 
of the phase noise is comparable to, or larger than, the variance of the random phase perturbation used 
in the one-bit feedback algorithm. This provides guidance on choice of system parameters such as LO 
quality, relaying delay, and cycle length. The open-source implementation of our prototype is publicly 
available, and hopefully provides a starting point for further implementation of solutions for distributed 
MIMO. 

There are many directions for future work. An important topic is generalization of our amplify- 
forward approach to provide scalable distribution reception over wideband dispersive channels. Possible 
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approaches include “filter-and-forward,” or amplify-forward on a per-subcarrier basis. Design challenges 
include timing synchronization and tracking schemes, and the development of parsimonious feedback 
strategies. Also, while our time division architecture yields implicit frequency synchronization, there 
may be many scenarios in which frequency division between long and short links is an attractive design 
choice, in which case explicit frequency synchronization is required. Finally, it is important to develop 
and evaluate designs that account for mobile nodes, possibly with different models addressing different 
potential applications. 


Appendix A 

Joint probability distribution of U and V 


in 13 


It was shown in 110] that for large N, the net effect of a random phase perturbation on the total signal 
can be modeled as shown in Figure |T4] The effect of phase noise (or random phase perturbations plus 



Fig. 14. Effect of a random phase perturbation on the total received signal 


phase noise) on the total signal can be modeled in an identical manner. Using equation (22) in 
large N and y, the following approximation can then be made: 



V6 « X5V + XfKj 
Vn « XnV + X<yi,n 
VSn w XSn.y 

The variables x<r } $, xy\. n and are zero-mean Gaussian random variables with variances tr^ & 


1 -X 2 5~P5K(y) 


2N ’ CJ 9t,n 2 N 

statistics of the random variables U = ys n ~ Vn and V = y$ - y are simply those of a bivariate Gaussian 
distribution, entirely characterized by the means of U and V, the variances of U and V, and the covariance 
between U and V. The development in the following subsections are similar to the development made 
in JfoIi -A ppendix C. 


l-Xn-PnK(y) 


and a kdn 


1 Xsn 2 ^ nK ^ , respectively. It can be seen that the joint 
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A. Mean and variance of U 

From Figure [l4l we can define the following terms 

1 N 

X5n = x<n t sn + jxj,6n = — E ( e j( - Si+n ^ - XSn ) 
- /V i =1 

1 N 

X n = Xfn,n + jX3,n = — E ^ ~ An) 

i =1 

1 N 

X6 = X<K,6+jX3,6 = — - Xs) 

iv i =1 

The mean and variance of t/ = y^ n - y n are then given by 

= ( XSn~Xn)y 

Var[t7] = E [(ax^ - x^n) 2 ] 


= ^E [(XSn + X* 5n - X n - X * n ) 2 

— — E + X^ n + X n + X n + ‘ 2 x§ n Xfi n + ‘ 2 x n X n — 2 x/j n X n — 2 Xfi n X n 


- 2 x* Sn x n - 2x* Sn x* n ] 


The first term of the previous equation is given by 

ekj= 4 eee[ 




XSn)(e- 


j(S,+n[) _ 


XSn) 


i=ll=l 


Since the phase perturbations Si and noise terms n\ have a symmetric distribution, it follows that the 
terms (e- 7 ^ i+n ^ ~X5n ) and ( e J+n >~X8n) have zero mean. Moreover, since the phase perturbations S t 
and noise terms n\ are independent for different nodes, the term in the sum are zero when i ± l and the 
previous becomes 

■ Nn] = jf- 2 EE [e^-(e^ + <> - X5n) 21 

1 N 

= -^Ee [^1 


E 


1=1 


where we used the approximation that, for small phase perturbations and/or phase noises, E ^ 
E [cos(<5j + n \)]. Similarly, one can compute 




E[x 5n x|J 


1 

iv 2 


N N 


EE 4 


(5f +n'i) _ 


XSn)(e 


~j( s i+ n 'i) _ 


X&nl 


i=ll=l 
N 


i JV r 

^S E 1 




X,5n ) 
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since E 


e j( s i+ n 'i) 


21 


X5n 


= 1 - Xs n (see 



). Following similar arguments, one can easily obtain the 


following expressions 


EK] = ~P„E E [^‘] 

1 V Z=1 

E K] = -^4 E[e ~^‘ ] 

E bn<] = -^(i-Xn) 


The crossed terms are computed as 


IE [Xfi n x n ] = ^ £ E [ 


N N 




XSn ) ( e - 


j(m) _ 


Xn) 


i=lZ=l 


The terms (e ? "( <5i+n «) ~xs n ) ar >d (e^ n ‘^ ~Xn) arc both zero-mean and independent variables, and therefore 
all terms in the summation equal zero, which leads to 


IE [. xs n x n ] = E [ x 5n x* n \ = E [x* Sn x n \ = E [x* 5n x* n \ = 0 


N 


The terms n(y) - jt T, Efe 2 ^^] depends on y only, and can be approximated by n(y) = e 4(1 for 


i =1 


>n y only. 

1 in [1QD; 


large y. The latter was derived in Conjecture 1 in 110]), where statistical mechanics arguments where used 
to determine that that the distribution of the phases around their mean follow an Exp-Cosine distribution. 


The variance of U can finally be written as 


Var[[7] = + 


2 N 


1 ~ Xn ~ PnK(y) 

2 N 


B. Mean and variance of V 

The statistics of V - ys - y can be deduced in a manner identical to the statistics of U, which leads 
to the following expressions: 


V] = (b-i)» 

Varm = 

L J 2 N 
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C. Covariance between U and V 

The two terms U and V are not independent, but the covariance between both can easily be obtained 
following similar arguments as before. 


CO \\U,V] = E[(.T<K i5n - X£H >n )xjH )5 ] 

= ^E [{x 5n + x* 5n - x n - x* ) (xs + X 5 )] 

= ^ E [x Sn XS + X &n x * 5 + X* 5n XS + X* Sn Xg ~ X n XS ~ X n X*g ~ X* n Xg - X* n X* 8 ] 

Note that the mean of the last four terms can immediately seen to be zero, since the phase perturbations 


terms and the noise terms are zero-mean and independent. The first term can be computed as follows: 

1 N N 

V[x8nX S ] 

iV j=ll=l 

- jp E E [«*'* (e j(26<+ni) - X8ne j5i - Xse j(Si+ni) + XSnXs)] 

1 N 

= Trt E E [ e2j<t>i ] (E[cOs(25j + Hi)] - XSnXS - XSXSn + XSnXS) 

^ i=l 

1 N 

= ^2 E E [ g2 ^] ( E [cOs(25 ?: + rii)] - XSnXS) 

Similarly, one can compute that 

1 N 

E [ X*Sn x S 1 = E E [e _2j01 ] (E[cOs(2<5j + m)] - XSnXS ) 

The crossed terms can be computed as follows 

1 N N 

E [xsnx* s \ = E E E - XSn)e~^{e~^ - xs) 

i=ll=l 

I N r 

i=l 

I N r 

= ^ E E - X«ne“^ + X<5nX5_ 

= ^ (Xn — XSnXs) 

and 

E [^n^l = T7 (Xn - XSnXs) 


Finally, the covariance between U and V can be written as 

C rjj tai _ Xn ~ XSnXS ~ P2SnXt(y) 

° VL ’ J ” 2N 


where p 2 g n = XSnXS ~ E[cos(2<5j + m)]. 
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